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Drug synergy allows a therapeutic effect to be achieved with lower doses of component drugs. Drug 
synergy can result when drugs target the products of genes that act in parallel pathways ('specific 
synergy') . Such cases of drug synergy should tend to correspond to synergistic genetic interaction 
between the corresponding target genes. Alternatively, 'promiscuous synergy' can arise when one 
drug non-specifically increases the effects of many other drugs, for example, by increased 
bioavailability. To assess the relative abundance of these drug synergy types, we examined 200 pairs 
of antifungal drugs in S. cerevisiae. We found 38 antifungal synergies, 37 of which were novel. While 
14 cases of drug synergy corresponded to genetic interaction, 92% of the synergies we discovered 
involved only six frequently synergistic drugs. Although promiscuity of four drugs can be explained 
under the bioavailability model, the promiscuity of Tacrolimus and Pentamidine was completely 
unexpected. While many drug synergies correspond to genetic interactions, the majority of drug 
synergies appear to result from non-specific promiscuous synergy. 
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Introduction 

Pathogens or tumors that do not respond to single treatments 
may be amenable to combined drug treatments (Yeh etal 2006; 
Azmi et al 2010} . Synergistic drug pairs have special potential 
for treatment, since they allow a desired effect to be achieved 
with a lower total dose of administered medicine (Greco et aU 
1995). Theoretical and experimental studies have shown that 
drugs that exhibit synergy for a specific effect are usually not 
synergistic for side effects (Lehar et aU 2009; Owens et aU 2010) . 
These considerations make synergistic drug pairs ideal 
candidates for treatment of pathogens or tumors. 

As a result of their medicinal potential, synergistic drug pairs 
are sought both experimentally and theoretically. Large-scale 
experimental drug synergy screens have found that synergistic 
drug pairs are rare (4-10%) (BorisyetaZ, 2003; Zhang etaZ, 2007; 
Farha and Brown, 2010). Therefore, methods which predict 
synergistic drug combinations have been developed, such as 
network analysis of expression profiles (Nelander etaU 2008), or 
clustering of chemogenomic profiles of drug perturbations 
(Jansen et al 2009). In addition, perturbations with multiple 



drugs have been used to understand the organization and 
robustness of biological systems (Lehar et al 2007) . 

One model for drug synergy is the parallel pathway 
inhibition model which suggests that two drugs will be 
synergistic if they inhibit two proteins on parallel pathways 
essential for an observed phenotype (Yeh et al, 2009) . To find 
protein pairs in parallel pathways, we used synergistic genetic 
interactions, where two combined genetic perturbations result 
in a phenotype that is more severe than expected (Tong et al, 
2004; Mani et al, 2008). Because a drug's action on a target 
gene product can yield an effect that is similar to a 
hypomorphic or loss-of-function allele of the gene encoding 
the drug target (Parsons et al, 2004; Hillenmeyer et al, 2008), 
two drugs may be synergistic for the growth inhibition 
phenotype if they target the products of two genes that have 
a synergistic genetic interaction. 

Several examples of genetically consistent synergy have 
been reported, such as the synergy involving a combination of 
either two mutations or two drugs that inhibit the DNA 
polymerase and thymidine kinase of Herpes Simplex Virus 
(Jacobson et al, 1995) . Also, drug synergy between Terbinafine 
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and Rapamycin is known to correspond to a synergistic genetic 
interaction between the target genes (Lehar et aU 2007} . 

Other models have been proposed to explain the mechanism 
of observed synergistic drug interactions (Jia et aU 2009} . The 
bioavailability model suggests that two drugs will be synergis- 
tic if one drug's action helps another drug's availability in the 
target cells (Zimmerman et al, 2007}, either by increasing the 
second drug's entry to the cell or by decreasing the second 
drug's degradation. For example, it has been reported that 1,8- 
Cineole enhances the effect of AZT, an anti-HIV drug, by 
increasing skin permeability (Narishetty and Panchagnula, 
2004} . Consistently, a number of drugs known to disrupt cell 
membranes, such as Triclosan in E. coli and ergosterol 
pathway-targeting drugs in yeast, have been found to exhibit 
many synergies (Ryder, 1992; Schweizer, 2001}. 

Combination therapies for which an underlying genetic 
explanation exists are attractive from a therapeutic standpoint. 
Such genetically determined drug synergies are expected to be 
more specific to particular pathways, thus having a reduced 
potential for unintended consequences. Concordantly, synergy 
which arises via the bioavailability model comes with the 
concern that bioavailability may be increased for many drugs 
beyond the one for which synergy is desirable, including drugs 
CO -administered for other indications. However, it may not 
always be possible to find an efficacious combination therapy for 
which the genetic underpinnings are known. Indeed, it is unclear 
whether correspondence between gene synergies and drug 



synergies, which we refer to as 'specific synergy', is the 
predominant form of drug synergy. 

If some drugs are inherently more likely to participate in 
synergistic combinations — for example, because they generally 
increase the bioavailability of other drugs — then the search for 
synergistic combinations might be accelerated by prioritizing 
combinations that contain promiscuous synergizers. Because 
'promiscuous synergy' by definition encompasses many different 
drugs, it may be relatively common when compared with specific 
synergy. 



Results 

We first established an experimental framework for classifying the 
interaction between two drugs as synergistic, independent, or 
antagonistic. For each drug pair examined, we examined all 
possible combinations of seven concentrations for each drug, such 
that the concentration for each drug ranged from 0 to its minimal 
inhibitory concentration (MIC} . Synergistic or antagonistic inter- 
action can then be classified according to significant deviation 
from non-interaction behavior under the Loewe additivity model 
of drug synergy (Greco et al 1995} . The Loewe additivity model is 
intuitively appealing, in that a drug does not interact with itself 
under this model. As viewed on a two-dimensional grid, a drug 
combined with itself will yield linear isophenotypic contours 
(isobols}, curves such that each point along the curve corresponds 




Growth 



Figure 1 Experimental set-up for classification of drug interactions. For each drug pair, growth rates were measured for all pairwise combinations of seven drug 
concentrations, linearly increasing from 0 to the minimal inhibitory concentration (MIC). 'Isophenotypic' curves describing drug concentration combinations which yield 
the same phenotype are shown on the left. Isophenotypic curves are expected to be parallel to the diagonal for independent drug pairs, concave for synergistic drug 
pairs, and convex for antagonistic drug pairs, according to Loewe additivity. For each drug pair, a drug interaction score (a) quantifying the concavity of the isophenotypic 
curve was computed, with a=0 defining independence and a taking negative or positive values when drugs are synergistic or antagonistic, respectively. Measurement 
error in a was assessed by examining the distribution of a for 25 self-self drug combinations, and drug pairs that had significantly smaller or larger a values were 
classified as synergistic or antagonistic, respectively. The white region shows the margin of variability between self-self interactions (-0.78 < a < 0.68, 95% confidence 
interval). Experimental examples of drug interactions are given on the right. In each example. Tacrolimus concentration is increased along the x axis. Drugs used 
in increasing concentrations along the y axis are Tacrolimus for self-self, Myriocin for independence, Latrunculin B for synergy, and Bromopyruvate for antagonism. 
For each example, sample isophenotypic curves are depicted in white. Yellow dots mark the longest isophenotypic curve, which was used to classify a drug interaction. 
See Supplementary Figure 1 for plots of all 200 painwise drug interaction tests conducted in the course of this study. 
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to a different drug concentration combination that yields the same 
growth phenotype. However, the contours will be convex when 
two drugs are synergistic, and concave when they are antagonistic 
(Figure 1). We used a mathematical model (see Materials and 
methods) to quantify the shape of isophenotypic contours where 
values near 0 indicate linearity (independence), negative values 
indicate convexity (synergy), and positive values indicate 
concavity (antagonism) . 

To measure experimental variation in this interaction score, 
we evaluated 25 'self-self drug combinations that would be 
expected to yield a score of 0 in the absence of measurement 
errors. This allowed us to classify drug interactions as 
independent, significantly synergistic, or significantly antag- 
onistic, depending on where the interaction score falls relative 
to the observed distribution of self-self interaction scores. We 
also considered Bliss and Gaddum's interaction models for such 
binary classifications (Berenbaum, 1989; Greco et al, 1995); 
however, we found that our model based on Loewe additivity 
provided the most intuitive and robust classification of drug 
interactions (see Materials and methods). Importantly, growth 
rates for self-self drug interaction experiments exhibited very high 
correlation along the symmetry axis (r=0.98, P<2.2 x 10"^^), 
indicating the reproducibility of our experimental system. 

To assess the overlap between synergistic genetic interac- 
tions and synergistic drug interactions, we focused on the 
yeast S. cerevisiae, the organism for which the largest number 
of synergistic genetic interactions is known. We compiled 



a catalog of 113 known chemical/target relationships in 
S. cerevisiae from the literature (Supplementary Table 1). All 
chemical compounds that yield a growth fitness defect when 
administered as single agents (referred to in this study 
as 'drugs') were considered, excepting compounds known 
not to act as inhibitors of their target. We next integrated the 
resulting drug-target catalog with known synergistic 
genetic interactions (Baryshnikova et al 2010; Costanzo et al, 
2010; Stark et al 2011; see Materials and methods). This 
yielded a set of 211 drug pairs predicted to be synergistic 
according to the parallel pathway inhibition model (Supple- 
mentary Table 2) . To assess the success of these predictions, 38 
of these drug pairs involving a total of 21 drugs were assessed 
for synergy. Except for three different drugs (DYC, FEN, and 
HAL) that were each used individually to inhibit Erg2, an 
ergosterol pathway protein, we used only one drug to inhibit 
each specific protein corresponding to a gene with a synergistic 
genetic interaction (see Table I for each drug's target name and 
function) . Nine of these thirty-eight tested drug pairs were also 
predicted to be synergistic according to the bioavailability 
model, as they included an ergosterol pathway-targeting drug. 

We found that 14 of the 38 tested drug pairs predicted to be 
synergistic using synergistic genetic interactions showed 
significant synergy (37% accuracy), while 11 of the 38 drug 
pairs (29%) showed significant antagonism (Figure 2). 
Although a set of negative control predictions tested within 
our framework would provide a better validation of the use of 



Table I List of all drugs used in this study, their abbreviation, highest dose used, target (bold if essential] and target function. Drugs that were found to be promiscuous 
are highlighted. Literature references are provided in Supplementary Table 1 



Drug 


Abbreviation 


Dose (|ig/ml) 


Target(s) 


Target function 


5 fluorouracil 


5FU 


28 


Cdc21 


Nucleotide metabolism 


Aureobasidin A 


ABA 


280 


Aurl 


Sphingolipid metabolism 


Amphotericin B 


AMB 


0.49 


Erg28 


Ergosterol metabolism 


Anisomycin 


ANI 


3.5 


? 




Benomyl 


BEN 


28 


Tubl,Tub2, Tub3 


Microtubule component 


Bromopyruvate 


BRO 


490 


Pykl 


Glycolysis 


CCCP 


C3P 


21 






Calyculin A 


CAL 


2.1 


Glc7 


Serine/threonine phosphatase 


Cantharidin 


CAN 


140 


? 




Chlorzoxazone 


CHL 


350 


? 




Cisplatin 


CIS 


80 


DNA 




Clozapine 


CLO 


105 


Coxl7 


Metallochaperone 


Cycloheximide 


CYC 


0.91 


Rpl41A, Rpl41B 


Ribosomal protein 


Dyclonine 


DYC 


49 


Erg2, Erg24 


Ergosterol metabolism 


Fenpropimorph 


FEN 


1.54 


Erg2, Erg24 


Ergosterol metabolism 


Haloperidol 


HAL 


56 


Erg2, Erg24 


Ergosterol metabolism 


Hygromycin 


HYG 


7 


Hocl 


Cell wall metabolism 


Latrunculin B 


LAT 


14 


Actl 


Actin 


Lithium 


LIT 


4500 


Inml, Met22 


Inositol metabolism, nucleotide 










metabolism 


Methotrexate 


MET 


1000 


Dfrl 


Nucleotide metabolism 


Methyl methanesulfonate 


MMS 


175 


DNA 




Myriocin 


MYR 


0.350 


Lcbl, Lcb2 


Sphingolipid metabolism 


Pentamidine 


PEN 


70 


Pntl 


Mitochondrial protein 


Quinine 


QNN 


1000 


Tat2 


Amino acid transport 


Quinomycin 


QMY 


30 






Radicicol 


RAD 


56 


Hsp82, Hsc82 


Chaperone 


Rapamycin 


RAP 


0.0049 


Tori, Tor2 


Protein kinase 


Staurosporine 


STA 


1.26 


Pkcl 


Serine/threonine kinase 


Tacrolimus 


TAC 


110 


Cnbl 


Phosphatase ^ 


Tamoxifen 


TAM 


2.8 


2 




Terbinafine 


TER 


10.5 


Ergl 


Ergosterol metabolism 


Tunicamycin 


TUN 


0.35 


Alg7 


Glycosylation 


Wortmannin 


WOR 


280 


Stt4 


Kinase 
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A 




STA RAD RAP STA RAP TUN RAD 




RAP WOR QNN LIT LAT CAL LAT 



] 

1 

Figure 2 Synergy tests between 38 drug pairs that target proteins encoded by synergistic genes. (A) Drug pairs were combined in 8 x 8 matrices where the 
concentration of one drug was linearly increased along each axis. The lowest concentration was 0 and the highest concentration was chosen close to the MIC for each 
drug. Fourteen drug pairs showed significant synergy (marked with S) while eleven were significantly antagonistic (marked with A). Thirteen drug pairs yielded scores that 
fell within the distribution observed for self-self drug pairs, and were classified as independent (unmarked). (B) Network representation of the drug interactions shown in 
(A). Edges reflect the interaction type between two drugs and the node pie charts represent the ratio of different types of interactions each drug has in this data set 
(green: synergy, white: independent; red: antagonism). For drug name abbreviations, see Table I. One example for drug synergy prediction using synergistic genetic 
interactions: MYR targets Lcb2 and QNN targets Tat2 is also shown. There is a known synergistic genetic interaction between LCB2an6 TAT2, which predicts that MYR 
and QNN will be synergistic according to the parallel pathway inhibition model of drug synergy. However, these drugs were found to be antagonistic. 




synergistic genetic interactions for predicting drug synergies, 
we initially noted that the prediction accuracy of the parallel 
pathway inhibition model was much higher than the 
previously reported background rate of antifungal synergy of 
4-10% (Borisy et al, 2003; Farha and Brown, 2010). 

Of particular note, all of the nine drug pairs that included an 
ergosterol pathway-targeting drug also exhibited synergy (100% 
accuracy). Indeed, these nine drug pairs — which represented 
only 24% of tested pairs — yielded 64% of the synergies among 
tested pairs, supporting the idea that drugs increasing bioavail- 
ability can provide an enhanced rate of synergy (P=1.2 x 10"^, 
Fisher's exact test). That ergosterol pathway-targeting drugs 
were enriched for synergy indicated that different drugs may 
have inherently different background synergy rates. In an effort 
to assess background synergy rates of different drugs, we tested 
all drug pairs in a square matrix of 13 drugs (Figure 3} . Of these 
13 drugs, 10 were partners in successfully predicted synergistic 
drug pairs in the initial screen (TER, FEN, HAL, DYC, PEN, TAG, 
LAT, BEN, STA, and RAP; see Table I for abbreviation key). Two 
drugs were partners in unsuccessful predictions (TUN and CAL) 
and BRO was added to the list arbitrarily. 

Of the 63 drug pairs in our test space that were not predicted to 
target synergistic gene products, we were surprised to observe 
that 21 (33%) showed synergy. The number of synergies each 
drug exhibited was not uniform (x^ test, P=0.024} , indicating the 
non-uniformity of background synergy rates for these 13 drugs. 
This observation suggests that there is no single universal 
background rate of synergy, which can be used as a basis for 
comparison of the success for a synergy prediction method. 

As an initial attempt to account for background rates of 
synergy, we considered the baseline synergy rate of each drug 
tested in our 13x13 matrix in turn. For each drug, we 
examined whether drug pairs involving that drug were more 
likely to exhibit synergy if they corresponded to a genetic 
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interaction. After correcting for multiple hypothesis testing, 
we found no evidence for significant enrichment of synergy 
among drug pairs corresponding to synergistic genetic inter- 
actions (using a threshold of 0.05 for adjusted P-values). This 
was despite the fact that the choice of drugs to fully test within 
the 13x13 matrix was initially biased to include drug pairs 
corresponding to genetic interaction. Thus, a consideration of 
baseline rates of synergy for each drug eliminates the 
significance of our observed correlation between synergistic 
genetic interaction and drug synergy. 

Where drug synergy does correspond to genetic interaction, 
we wished to further analyze whether the observed drug 
synergies are the result of synergistic genetic interactions 
between the drug targets. Therefore, we conducted drug 
interaction experiments with an S. cerevisiae strain for 
which the CNBl gene — encoding the target of TAG (Parsons 
et aU 2004}— was deleted. CNBl has a synergistic genetic 
interaction with ACTl and PKCl (targets of STA and LAT, 
respectively), and we have shown that drug pairs TAG-LAT 
and TAG-STA are highly synergistic. We hypothesized that if 
the synergies between these drugs arise from the same 
mechanism as the synergistic genetic interactions, then these 
drug synergies should not exist in a cnbl null mutant. 
However, we observed both of these drug synergies in cnhlA 
strains. This suggests that the observed synergies of drug pairs 
TAG-LAT and TAG-STA, successfully predicted on the basis of 
genetic interaction, are not mediated solely by the CNBl- ACTl 
and CNBl-PKCl synergistic genetic interactions. The observa- 
tion that TAG exhibits synergies even when its known target 
has been deleted also suggests that there are other targets of 
TAG in S. cerevisiae. These results exemplify the idea that the 
existence of multiple targets for each drug can comphcate tests 
of the correspondence between genetic interaction and drug 
synergy. 
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TER PEN TAG FEN HAL DYC LAT STA BEN CAL RAP TUN BRO 




Figure 3 Synergy tests between all pairs among 1 3 drugs. (A) Drug pairs were combined in 8 x 8 matrices where the concentration of one drug was linearly increased 
along each axis. Green squares correspond to drug pairs that target proteins encoded by synergistic genes. We found 32 significant synergies (S) and 27 significant 
antagonisms (A) among the 78 drug-drug interactions tested. Hierarchical clustering of drug interaction score profiles is shown on the left. (B) Network representation of 
the synergistic and antagonistic drug interactions shown in (A). Edges reflect the interaction type between two drugs and the node pie charts represent the ratio of 
different types of interactions each drug has in this data set (green: synergy, white: independent; red: antagonism; independence edges are omitted for clarity). Gray box 
indicates the promiscuous synergizers learned from the drug interaction network in this figure. Yellow circle shows the drugs that target the ergosterol pathway in yeast. 
For drug name abbreviations, see Table I. 



The observation that drugs have different background 
synergy rates suggested a strategy whereby pairs that involve 
at least one of these promiscuous synergizers are prioritized for 
testing. The number of synergies shown by each drug is a 
proxy for their background synergy rate. However, each of 
these synergies may be due to the high background synergy 
rates of the partner drugs. For a more robust analysis of 
background synergy rates of drugs, we clustered the interac- 
tion score profiles of the drugs in the 13 xl3 matrix. Not 
surprisingly, the three Erg2 -targeting drugs were clustered 
together. Interestingly, these three and the other ergosterol 
pathway-targeting drug TER were clustered with TAG and 
PEN, despite the fact that the latter two drugs are not known to 
disrupt the cell membrane. All pairwise interactions within 
this set of six drugs were synergistic, except between Erg2- 
targeting drugs (which may be explained by competition 
between these drugs at the same target site). We refer to 
members of this cluster as promiscuous synergizers. 

We also carried out a network analysis using the drug 
interaction network shown in Figure 3B in order to determine 
the smallest set of drugs that could 'explain' all observed drug 
synergies. When two drugs acted synergistically, we attributed 
the drug synergy to the drug with the highest level of background 
synergy. By individually removing the drugs with the highest 
background synergy rate and recomputing the number of 
synergies each drug exhibits in the resulting network, we found 
that the parsimonious classification of six drugs as promiscuous 
explained all the drug synergies in this data set. This smallest set 
coincided exactly with the promiscuous synergizers found by 
clustering analysis. Because the remaining seven drugs in our 
analysis exhibited relatively few synergies, we refer to these as 
chaste synergizers. 

© 201 1 EMBO and Macmillan Publishers Limited 



We next wondered whether promiscuity and chastity of 
drugs are intrinsic drug properties that would predict drug 
interaction behavior outside of the 13 xl3 drug matrix in 
which these properties were learned. Therefore, we further 
examined three promiscuous (TER, PEN, and TAG) and three 
chaste drugs (LAT, STA, and BEN), assessing synergy for each 
of these six against an array of 14 arbitrarily chosen drugs 
(Figure 4). Fourteen of sixteen synergies observed in this 
screen included a promiscuous synergizer drug, indicating a 
significant enrichment (Fisher's exact test, P=8.7 x 10""^} and 
supporting the idea that promiscuity is an intrinsic property with 
predictive value. Remarkably, LAT and STA did not show synergy 
with any of the 14 drugs, supporting the idea that they are 
intrinsically chaste. We found that BEN showed synergy with 
two arbitrarily chosen drugs (G3P and QMY) and it clustered 
with the other three promiscuous drugs in this new data set. 
While further experiments will be required to confidently 
characterize BEN as either promiscuous or chaste, these 
observations support the inherence of promiscuity for TER, 
PEN, and TAG. These results further bolster the conclusion that 
background synergy rates of different drugs are highly non- 
uniform. 

We next wondered whether drugs have intrinsic tendencies 
towards antagonism, as well as synergy. An analysis of the 
data in the 13x13 square matrix of drug interactions does not 
reject the null hypothesis that background antagonism rates of 
drugs are uniform (x^ test, P=0.097). However, we note that 
BRO, which showed antagonism with every drug against which 
it was tested, was significantly enriched for antagonisms 
iP=2.9 X 10"^}. Moreover, the synergy and antagonism rates of 
individual drugs have a significant negative correlation 
(r=-0.78, P=1.6xl0"^}. Thus, the antagonism background 
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Figure 4 Synergy tests between 6 drugs (3 predicted to be promiscuous synergizers and 3 predicted to be chaste synergizers) and 14 arbitrarily chosen drugs. (A) 
Drug pairs were combined in 8 x 8 matrices where the concentration of one drug was linearly increased along each axis. Promiscuous synergizers (top three rows) 
showed 14 synergies (S) while chaste synergizers (bottom three rows) showed 2 synergies. Thirty-seven significant drug antagonisms (A) were observed, predominantly 
involving chaste synergizers. Green squares correspond to drug pairs that target proteins encoded by synergistic genes. Hierarchical clusterings of drug interaction score 
profiles are shown on the left and top. (B) Network representation of the synergy and antagonism drug interactions shown in (A). Edges reflect the interaction type 
between two drugs and the node pie charts represent the ratio of different types of interactions each drug has in this data set (green: synergy, white: independent; red: 
antagonism; independence edges are omitted for clarity). Gray circle indicates the promiscuous synergizers learned from the drug interaction network in Figure 3B. For 
drug name abbreviations, see Table I. 



rates of individual drugs may be a useful exclusion criterion for 
accelerating the search for synergistic drug combinations. 



Discussion 

In this study, we systematically explored drug combinations to 
evaluate specific and promiscuous models of drug synergy. We 



developed a sensitive drug interaction classification assay and 
used it to find 38 synergistic drug pairs, of which only 1 
(Terbinafine-Rapamycin) had been previously reported 
(Lehar et aU 2007). 

To identify cases of specific synergy, predicted by synergistic 
genetic interaction, we integrated a newly assembled catalog 
of drug targets with the synergistic genetic interaction 
network. Although we tested 38 drug pairs corresponding to 
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genetic interactions and found 14 (37%) to be synergistic, it is 
difficult to assess whether these drug synergies are better 
explained by genetic interaction or promiscuous synergy. 
Indeed, after accounting for background rates of synergy for 
each drug, our results did not show a significant enrichment 
for drug synergy among drugs targeting the products of 
genetically interacting genes. 

There are many potential explanations for this observation. 
First, either or both of the drugs may not inhibit their targets 
completely or at all. Although chemogenomic studies have 
shown that the deletion of a target gene is often similar to 
treatment by the drug (Hillenmeyer etal 2008), and although we 
excluded drugs that were known not to act as inhibitors, not all 
of the drugs within our tested set are known to act as inhibitors. 

Second, each drug may have multiple targets. Action on 
additional targets may raise the background rate of drug 
synergy, making it more difficult to observe a significant 
enrichment for synergy predicted via genetic interaction. In 
rare cases, action on these additional targets may even 
suppress the effects of action on the known targets used to 
make drug synergy predictions. 

Third, some drug targets in our literature-derived catalog 
may have been erroneously identified. The misidentification of 
targets would naturally have a profound effect on the apparent 
overlap between genetic interaction and drug synergy. 

Fourth, intrinsic promiscuity of individual drugs may be the 
dominant factor in predicting drug synergy, again making it 
difficult to discern significant enrichment for synergy for drug 
pairs that correspond to genetic interactions. Indeed, our 
results support the idea that mechanisms other than parallel 
pathway inhibition of drug synergy are largely responsible for 
the synergistic drug interactions that are observed in our study. 

Another model of drug synergy proposes that drugs may 
exhibit synergy if a drug's action helps a second drug's 
availability in the cell [hioavailahility model). This model was 
systematically evaluated in our study, by testing synergy of drug 
pairs that included drugs that are known to destabilize the cell 
membrane by targeting the ergosterol pathway in yeast. We have 
shown that such drug pairs are enriched for drug synergy. We 
further classified drugs that partake in numerous synergies as 
'promiscuous synergizers' (Farha and Brown, 2010). Interest- 
ingly, two of the drugs we classified as promiscuous synergizers 
are not known to destabilize the membrane. Further studies will 
be required to determine whether there are other classes of 
promiscuous synergizers that do not correspond to the bioavail- 
ability model of drug synergy. In fact, after > 50 years since the 
bioavailability model for drug synergy was suggested (Jawetz 
and Gunnison, 1953), direct experimental evidence for this 
mechanism has not been produced. However, the analyses 
presented here strongly suggest that promiscuous synergizers 
enhance the potency of other drugs. Therefore, promiscuous 
drugs have a potentially valuable use in accelerating the search 
for synergistic drug combinations. 

Additional models explaining drug synergy have been 
proposed. The physical interaction model proposes that two 
drugs physically interact to make a more potent compound. This 
model, which is largely discredited (Jawetz and Gunnison, 
1953), was not evaluated here. The same target model proposes 
that two drugs will be synergistic if they target different sites on 
the same protein (Krogstad and Moellering, 1986). Incidentally, 
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it has also been proposed that two drugs will be antagonistic if 
they target similar sites on a protein (competitive binding model; 
Krogstad and Moellering, 1986). In this study, we tested all 
pairwise combinations of three drugs that target the same 
protein, Erg2. Two of these drug pairs were found to be 
antagonistic, while none of them exhibited synergy. These 
observations support the idea that these three Erg2 binding 
drugs, as predicted by the similarity of their chemical structures, 
bind to similar pockets on the target protein (Giaever etal 2004} . 

Our results show clearly that background synergy rate is 
drug dependent. This observation is critical for the assessment 
of any drug synergy prediction method, since it should account 
for the background synergy rates of all drugs that are in a drug 
combination. 

The wealth of drug interaction results that are obtained in 
our study allows us to comment on mechanisms of drug 
antagonism. Two models for drug antagonism have been 
proposed. The first, referred to as the competitive binding 
model, was mentioned above with the same target model for 
drug synergy (Krogstad and Moellering, 1986). The other drug 
antagonism model suggests that drug antagonism occurs 
when a drug's action impedes the metabohsm of the cell so 
that the second drug cannot exhibit its optimal effect 
[metabolic imbalance model; Bollenbach et aU 2009; Xavier 
and Sander, 2010} . Concordantly, we note that Bromopyruvate, 
which is known to inhibit glycolysis, was antagonistic with all 
drugs tested, including promiscuous synergizers. Thus, our 
results hint that background drug antagonism rates are also 
non-uniform. Notably, a previous study suggested, in agree- 
ment with the metabolic imbalance model for drug antagon- 
ism, that drugs that target unrelated pathways are more likely 
to show antagonism (Lehar et aU 2007} . 

The study presented here provides a systematic exploration 
of drug synergies in the context of mechanistic models of drug 
synergy. We found 14 drug pairs to be consistent with specific 
synergy, that is, for which genetic interaction provides some 
evidence that the drugs target parallel pathways. Our results 
also suggest that discovery of synergistic drug pairs can be 
accelerated by considering the background synergy rate 
individual drugs. Although the predominant form of drug 
synergy in our study is promiscuous synergy, it remains 
possible that genetic interaction will be a good predictor of 
drug synergy in other contexts. For example, genetic interac- 
tion might predict synergies between chaste synergizers. 
Proper study of this question awaits an expansion in the 
number of known drug/target relationships and a further- 
expanded map of genetic interactions. 

Materials and methods 
Experimental materials and methods 

5FU, ANI, BEN, BRO, C3P, CHL, CIS, DYC, FEN, HAL, LIT, MET, MMS, 
PEN, QNN, TAM, and TER were purchased from Sigma- Aldrich. AMB, 
CAL, CAN, CLO, CYC, HYG, LAT, MYR, QMY, RAD, RAP, STA, TAC, 
TUN, and WOR were purchased from AG Scientific. See Table I for 
abbreviation key. ABA was purchased from Clontech. All drugs were 
dissolved in DMSO (except Cisplatin which was dissolved in DMF) and 
kept at -20°C (except Cisplatin and Bromopyruvate which was freshly 
prepared due to loss of activity upon freezing). All experiments were 
conducted with S. cerevisiae strain BY4741, purchased from Open 
Biosystems. Yeast cells were grown in YPD (1% yeast extract, 2% 

Molecular Systems Biology 201 1 7 



Systematic exploration of synergistic drug pairs 
iVI Col^ol et al 



bacto-peptone, 2% glucose) overnight (ODeoo^S) and diluted to an 
ODeoo of 0.01 in YPD with the desired drug concentrations controlled 
for solvent concentrations. Cells were grown for 24 h in 96-well plates 
in a Tecan F200 microplate reader, with OD595 readings of cell density 
recorded every 15 min. For each drug-drug interaction assay, we took 
6144 cell-density measurements (detailed growth time courses for each 
of 64 concentration combinations with 96 time points in each time 
course). Synergy assays were performed for 200 drug pairs in total, so 
that a total of > 1 million measurements were performed during this 
study. We used the area under growth curve of each condition as a 
metric of cell growth, after discarding the first 10 measurement points. 
The entire collection of growth measurement data is available as 
Supplementary Data Set 1 . 



Synergy assessment 

To assess synergy, we used the Loewe Additivity Model, which defines 
a drug's interaction with itself as no interaction (for a comparison with 
other drug synergy models, see below) (Greco et al, 1995). Eight 
different concentrations of two drugs were combined in a two- 
dimensional grid, where the lowest concentration for each drug was 0 
and the highest concentration for each drug was close to the MIC. The 
concentration of one drug was linearly increased on each axis (see 
Table I for the highest dose used for each drug) . 

The marginal effect of each drug alone was determined by assessing 
the growth under 0 concentration of the other drug. Additional grid 
points were imputed (using Akima interpolation; Akima, 1970) and 
contours were plotted, connecting the marginal growth rates of each 
drug through the grid along 'isophenotypic' curves. The maximum of 
the lowest marginal growths of each drug was used as a threshold, 
where any growth lower than this value was not considered (as 
contours could no longer connect the two marginal axes) . The grids 
were then rescaled from 0 (low growth) to 1 (highest growth) and the 
longest isophenotypic curve (i.e., the curve described by the greatest 
number of data points) was selected for assessment. The following 
model was used to summarize this curve: 




where x and y are the normalized drug concentrations and a is a free 
parameter describing the 'bend' of the curve. When a=0, the contour 
connecting the two marginal growth rates is linear (e.g., describing the 
case when a drug is tested against itself). When a<0 or a>0, drug 
synergy or antagonism is concluded, respectively. 

To assess measurement error, 25 drugs were experimentally tested 
for interaction with itself (since by definition, there is no synergy for a 
drug with itself, any deviations of a from 0 represent measurement 
errors) . We modeled the a values of drug self-interaction experiments 
as normal (as we did not find sufficient evidence to believe otherwise, 
KS test, D=0.15, P=0.58), with a mean of fi of -0.05 and standard 
deviation a of 0.37. Hence, two drugs were considered synergistic or 
antagonistic when their a value was significantly smaller (a < -0.78) 
or larger (a > 0.68) than the experimental mean of the drug self 
interaction scores, respectively. 

The synergy levels between drugs were also assessed using Bliss 
Independence and Gaddum's non-interaction models. For Bliss 
Independence, parameter p that minimizes the following metric was 
computed to assess the overall synergy between drugs 1 and 2: 

^[f(drugl[^] +drug2[3;]) - pxf(drug l[^])xf(drug2[3;])]^ 

and for Gaddum's non-interaction model, parameter y that minimizes 
the following metric was computed: 

[f(drug 1m + drug2[y]) - yx max{ f( drug 1 f(drug2[3;])}]^ 

where f(drug 1 [x]] represents the cell growth in drug 1 at concentration 
X, and p and y are free parameters describing the deviation from null 
model surfaces. According to these drug interaction models, p=l or 
y=l represents non-interaction between two drugs. When these 
parameters are smaller or larger than 1, synergy or antagonism is 
concluded, respectively. Parameters p and y for self-self drug 
interactions were approximately normally distributed. A comparison 
of the estimated p values with a scores estimated above showed that 



these two parameters had a very high correlation (r=0.79, P< 10"^^); 
however, y values showed correlation with neither a nor p parameters. 
Moreover, according to the estimated y values for our experiments, 
none of the tested drug pairs showed synergy or antagonism. While a 
parameters found only 1 self-self interaction as interacting (P<0.05), 
the Bliss Model found 2 self-self interactions as interacting (P>0.05). 
In the light of these observations, we concluded that Bliss Indepen- 
dence and Gaddum's non-interaction models, while useful for finding 
synergies between specific concentrations of two drugs, does not give 
robust results for an absolute measure of synergy between two drugs. 
Hence, we used the Loewe additivity model as our drug synergy- 
antagonism classifier. 



Genetic interaction data 

In addition to published sources of genetic interaction data (Baryshnikova 
et al, 2010; Costanzo et al, 2010; and Stark et al, 2011) we also used 11 
unpublished interactions from the Boone and Andrews Labs: AURl- 
MET22, ERG2-ACT1, H0C1-T0R2, ACTl-TORl, ACTl-CNBl, INM1-T0R2, 
INMl-CNBl, LCB2-TAT2, LCB2-HSP82, LCB2-CNB1 HSC82-PKC1. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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